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Abstract 

The establishment of protected areas or parks has become an important tool for 
wildlife conservation. However, frequent occurrences of human-wildlife conflict at the 
edges of these parks can undermine their conservation goals. Many African protected 
areas have experienced concurrent declines of apex predators alongside increases in both 
baboon abundance and the density of humans living near the park boundary. Baboons 
then take excursions outside of the park to raid crops for food, conflicting with the human 
population. We model the interactions of mesopredators (baboons), apex predators 
and shared prey in the park to analyze how four components affect the proportion of 
time that mesopredators choose to crop-raid: 1) the presence of apex predators; 2) 
nutritional quality of the crops; 3) mesopredator “shyness” about leaving the park; and 
4) human hunting of mesopredators. We predict that the presence of apex predators 
in the park is the most effective method for controlling mesopredator abundance, and 
hence significantly reduces their impact on crops. Human hunting of mesopredators is 
less effective as it only occurs during crop-raiding excursions. Furthermore, making crops 
less attractive, for instance by planting crops further from the park boundary or farming 
less nutritional crops, can reduce the amount of time mesopredators crop-raid. 


Keywords: human-wildlife conflict, mesopredator shyness, crop subsidies, crop-raiding, apex 
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1 Introduction 


In many African protected areas (hereafter “parks” 
tors, mesopredators, and prey has been shifting, 
leopard, spotted hyena) are disappearing due to poaching (Brashares et al. ^ 2004 Kaltenborn 


the balance between humans, apex preda- 
In most parks apex predators (e.g. lion. 


et al., 2005), disease (Murray et al. ^ 1999), and habitat conversion or reduction (Balme et al. 


2010). Simultaneously, the density of human populations around parks is increasing (Joppa 


et al, 2009 Wittemyer et al, 2008), compounding pressure on apex predators and increasing 


mortality due to poaching and human-wildlife conflict (Lindsey et al, 2005 Woodroffe et al 


2007). In many parks their prey, such as ungulates, are also subject to high rates of hunting for 
meat (Thirgood et al, 2004). While apex predators and their prey are decreasing, in contrast 


populations of mesopredators, such as baboons, are increasing. For example, over the period 
1968 - 2004 all large apex predators became extinct in 3 of the 6 parks studied by the Ghana 
Wildlife Division, but baboons had a 365% increase in observations and 500% increase in range 
(Brashares et al, 2010). This increased abundance of baboons, which occurs in many parks 
throughout Africa, frequently results in crop-raiding as the populations spill over into farmed 
land outside the park (Hill, 1997, 2000). Baboons are able to exploit all trophic niches. They 


can act as predators, compete for browse with ungulates and livestock (Strum & Western 


1982), and exploit domesticated landscapes. It has been shown that baboons respond quickly 


to newly available resources in terms of fecundity (Bercovitch & Strum, 1993), and have a high 
rate of potential demographic increase for a large primate. Given sufficient resources baboon 
populations could potentially increase at roughly 20% annually (see Appendix A). 

Increased baboon abundance can pose a serious problem for people living at or near the 
boundaries of parks because of the potential for a large percentage and wide variety of crops 


1997 

1998 

Tweheyo et al 

2005 

). This destruction of crops is detrimental to human livelihood 


and education (due to children and women staying in the helds to defend the crops (Mackenzie 


et al , 2015)). It can also contribute to negative attitudes toward the park, potentially under¬ 


mining conservation aims (Hartter et al, 2014, Ryan et al , 2015; Tweheyo et al, 2005). Thus, 
understanding why baboons crop-raid, and, importantly, how to control levels of crop-raiding, 
could be very benehcial both to park conservation goals and the livelihoods of people living 
near them. 

In a recent study, Nishijima et al (2014) presented a mathematical model of mesopredator 
release when the mesopredator has both prey shared by an apex predator and alternative, 
unshared prey. In particular, they examined the effect of including alternative prey on the 
change in shared prey abundance when apex predators are lost. Although mesopredator release 
does not necessarily have a negative effect on the shared prey species ( [Brashares et al , 2010), 
a large supply of alternative prey can intensify mesopredator release and its deleterious effects 
on the ecosystem. 

We expand on the framework developed by Nishijima et al (2014) to understand the 
influence of human subsidies on this crop-raiding mesopredator dynamic. In this case the 
subsidy is crops, which exist outside the boundaries of parks. The continual availability of 
crops could enhance baboon population increases and thus feed back into more crop-raiding. 
However, unlike the unlimited alternative resource of Nishijima et al (2014), in this model 
we assume that the boundaries of the park pose an added consideration for baboons. The 
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presence of humans, domestic animals, and reduced cover may make baboons reticent about 
leaving the park, even while increased baboon abundance within the park increases pressure 
to leave the park to crop-raid. By incorporating a more complete mathematical description 
of this behavioral trade-off in potential baboon crop-raiding, we can link previous work on 
mesopredator release with human subsidy control and management strategies, to gain an 
understanding of the balance in managing the damaging effects of crop-raiding. 

In African parks, many animals crop-raid including multiple primate species (e.g. gorillas 


(Hill 


(Biryahwaho, 2002 Hill, 1997); chimpanzees (Madden 1999); redtail, colobus, vervet, blue 


1997) and golden (Biryahwaho, 2002[ ) monkeys), ungulates such as buffalo, bushbuck, 
and duikers (e.g. Plumptre et al. (1997)), and elephants (e.g. Naughton-Treves & Treves 
(2005); Osborn (2004)). These animals may also face behavioral barriers to leaving the park to 
crop-raid, and will in many cases be affected by predation and/or human hunting. Therefore, 
our analysis can extend to a broad range of crop-raiding species rather than focused purely 
on baboons. Additionally, mesopredator release with human subsidies is not limited to parks 
landscapes in Africa. For example, dingo abundance in Australia is carefully controlled due 
to their predation upon farmed cattle. However, this has led to mesopredator release of foxes 
and cats and concurrent declines of their prey; foxes often beneht from human subsidies as 


sheep farming leads to an abundance of rabbits (Johnson et al.. 2007). Whilst this is not a 


direct analogue to our system, parallels can be drawn to give insight into a wider spectrum of 
ecosystems. We thus present a fairly specihc case to illustrate a framework that is flexible to 
many scenarios. 

In this paper we use a mathematical model to explore how predation, human hunting, 
mesopredator shyness, and quality of crops can all interact, leading to different proportions of 
time spent crop-raiding by a mesopredator. Model parameters are chosen with an eye towards 
the particular case of baboons, although the model is more general. Our model could be 
adapted to any scenario in which there are two competing species with a shared predator, 
where one of the competitors has access to a separate source of food such as human subisides 


(Oro et al., 2013). We explore the sensitivity of the dynamics to the parameters (including 


the willingness of the baboons to leave the park), and to changes in three potential controls: 
hunting of apex predators, hunting of mesopredators, and quality of the food subsidies. We use 
this analysis to assess how crop-raiding may be reduced by human or apex predator control. 


2 Model 

2.1 Model Outline 

We model the population abundances of mesopredators (baboons), apex predators, and shared 
prey (ungulates) as they interact with each other, their shared resources, and with human 
hunting. We use the term ungulates to represent many potential species of ungulates within 
the park. A schematic for the modeled ecosystem is provided in Figure [T] This shows that apex 
predators, at the top of the food-chain, prey upon both the baboons and the ungulates. The 
park itself provides resources for the baboons and ungulates, both in terms of food availability 
and living space. We do not model this resource explicitly, but rather implicitly through the 
carrying capacity of baboons and ungulates, i.e. how many of these animals the park is able 
to sustain. The baboons and ungulates compete for this shared resource; their impact on each 
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Figure 1: The interactions between the apex predators, baboons and ungulates with their 
available resources and the role of human hunting. Variables directly included in the model 
are shown as blue rectangles, orange ovals represent resources which impact carrying capacities, 
and the green hexagon shows the impact of human hunting. Baboons spend a proportion of 
their time, t{B) living off the park resources and the rest of their time, 1 — crop-raiding. 
Baboons are only hunted by humans while crop-raiding, and baboons and ungulates share an 
in-park resource pool, which determines carrying capacity. 


other arises from how well they exploit the shared resource. 

Baboons have an additional resource of agricultural crops, inaccessible to ungulates, which 
increases the number of baboons that can be sustained within this system as they are only 
restricted by spatial and social requirements (Warren et al, 2011). Furthermore, agricultural 
crops are often high in nutritional value therefore baboons are able to reproduce at a higher 
frequency when they feed off crops compared to the park resources. 

The human population outside the park is also implicit. We need only consider its impact 
on the apex predators, baboons, and ungulates. In this model humans hunt all three animal 
groups. However, they only hunt baboons because they raid agricultural crops, thus hunting 
only occurs when the baboons are crop-raiding. 

The model assumes that baboons spend a portion of their time, t{B), living off the re¬ 
sources of the park, in competition with ungulates, and the rest of their time crop-raiding, 
1 — t{B). We assume that crops are available for baboons through the entire year. We trans¬ 
late this into a mathematical model (Eqns. ([^-(|^) where B{t), S{t) and P{t) are the numbers 
of baboons, ungulates and apex predators at time t respectively: 


4 









dt 

d^ 

dt 

dP 

dt 
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=^BdpfBp{B)P + rfepil — 6p)fsp{S)P — (/i + Hp)P. 


( 1 ) 

( 2 ) 

(3) 


Both baboon and ungulate populations experience logistic growth with growth rates {rB,rs) 
and carrying capacities {Kb, Ks) determined by the available resources in the park. However, 
this is regulated by competition, with aBs being the rate of competition by S on B. The 
baboons only spend a proportion of their time t{B) feeding on park resources and competition 
between baboons and ungulates only occurs during this time. Hence, the competition term 
in (|^ is multiplied by t{B). When baboons are feeding on agricultural crops the population 
growth is still determined by a logistic term, but their growth rate is increased by rate a and 
the carrying capacity by rate fd to represent the benehts of crop-raiding - enhanced nutritional 
quality of food and reduced limitations on availability. All three populations are affected by 
human hunting {Hb, Hs, Hp), but baboons only when they are crop-raiding, hence Hp is 
multiplied by 1 — t{B) in Q. Apex predators survive and reproduce by feeding upon both 
baboons and ungulates, spending a proportion dp oi their time hunting baboons and the 
rest of their time hunting ungulates. Apex predators are also free to roam outside of the 
park and hence they are able to hunt baboons at all times regardless of whether baboons are 
crop-raiding (the predation terms do not contain t{B)). For now we use a general functional 
form {fBp{B),fsp{S)) to represent the predation terms, which could change for different 
ecosystems. For a full description of the parameters and their values, see Table [T] 

There is some evidence that baboons prey upon ungulates opportunistically, such as Thom¬ 
son’s gazelles in Gilgil, Kenya (Bercovitch & Strum, [1993 ; Strum & Western, 1982) and on 
Kob and goats, in and near parks in Ghana (J. Brashares Pers. Obs.). This conforms to a 
mesopredator release scenario in that loss of the top apex predators allows baboons to become 
the mesopredator, controlling ungulate population abundance both through competition for 
resources and by predation. We model this indirectly by imposing a large competitive effect of 
baboons on ungulates. In our model the ungulate population represents a number of ungulate 
species, as well as potentially smaller primates, all of whom are in competition with the ba¬ 
boons for park resources. Since baboons will not prey upon all species of ungulates we do not 
model direct baboon predation upon this category but rather subsume it into competition. 


2.2 Time Spent Crop-Raiding 

To understand the proportion of time that baboons focus on park resources, t{B), we have to 
make some assumptions about the factors determining when baboons crop-raid. We assume 
baboons choose optimally how to split their time between crop-raiding and park resources, 
based upon the benehts and costs of each action. Built into this is the idea that baboons will 
always spend some portion of their time feeding off park resources. The baboon population 
can grow with rate rp when utilizing park resources and at rate rpa when crop-raiding. 
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Parameter Interpretation Value 
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parameter values can be found in §2.3 












However, they are susceptible to human hunting when crop-raiding at rate Hb- If we did not 
include any impact of shyness, previous studies on patch choice ( Krivanf 1997) indicate that 
baboons would optimally split their time to feed proportionally upon park resources and upon 
crops in the ratio —— 5 — : rga-Hg respectively. However, we additionally consider that 
their level of crop-raiding is based upon a crowding parameter 7 (see Table baboons are 
”shy” about leaving the park, but higher baboon numbers encourages them to do so. This is 
introduced by considering the number of baboons in relation to their carrying capacity and 
how shy they are, i.e. ^ — 7 . This leads to the following form for t{B): 


t{B) = 


Tb 


tb + TBa - H 


+ 


B 


TBa - 


B 


Tb + TBa - 



l + exp(cr(£- 7 )) 


(4) 


The exponential term allows a smooth transition from 100% park resources to a combination 
of resources and crop-raiding, as B increases. The hrst term in this equation is the minimum 
proportion of time that baboons will spend feeding upon park resources, which is supple¬ 
mented by the second term when either shyness is high or baboon numbers are low. As B 
increases, the last term will approach 0 which indicates that baboons are spending the mini¬ 
mum proportion of time only in the park. Hence, the proportion of time crop-raiding would 
approach rga-HB _ 

rg+rgOl-HB 

Similarly, we assume that apex predators split their time predating upon baboons and 
ungulates based upon the relative benehts from each source, which leads to 


6p = 


^BfBp{B) 


fBp{B) 


^pfppiB) + rjeBfspiS) fBpiB) + rifsp{S) 


(5) 


for the proportion of time the apex predators spend predating baboons (and 1 — for ungu¬ 
lates). eB and are the rates determining how efficient the apex predators are at converting 
prey (baboons and ungulates, respectively) into reproduction. The functional forms for the 
predation of baboons and ungulates by apex predators are taken to be of Holling Type H 
form, such that 

^ ra’ 

where rriB is the maximum predation rate on baboons and is the number of baboons at 
which predation is half its maximum. A similar equation holds for fspiS) (see Table [^. 


2.3 Parameter Estimates 

Parameters are estimated from literature where possible, with the parameters calibrated to 
produce dynamically realistic scenarios for co-existence of baboons, ungulates and apex preda¬ 
tors (Table [^. We use available data from Greater Addo Elephant National Park, South 
Africa, which includes forested parkland but scale it to a smaller, more average park size. 
Since we assume that the S(t) class can represent many species of ungulates, and even small 
primates, with which the baboons compete for park resources, we add together estimates of 
ungulate abundance for the carrying capacity and use average birth and death rates from 
multiple species to calculate the growth rate. Growth rates for baboons are found in |Smuts~^ 


Nicolson (1989) and Ryan (2015) and carrying capacities in Hayward et al. (2007), with the 
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higher estimates used to quantify the increases due to crop-raiding. Growth rates for ungulates 
are found in Spinage (1972) and carrying capacities in Hayward et al. (2007). Estimates for 
apex predator longevity are found in Paemelaere & Dobson (2011) as well as information on 
fecundity rates to inform parameters and p. Estimates for the maximum predation rates 
{rriB, ms) and half-saturation constants {Gb, Gs) were taken from Altman & Dittmer (1968). 


Average population size of apex predators within different parks is found in Hayward et al. 


(2007), which was used to confirm the population estimates were reasonably accurate. The 
parameters chosen are intended for a conceptual understanding of how the different species in¬ 
teractions affect baboon crop-raiding generally. A more in-depth analysis at a specific location 
would require species data at that location. 


2.4 Method Outline 

We are interested in the effect of varying the parameters on the amount of time the meso- 
predator is expected to crop-raid. We initially take a local stability approach, whereby we 
numerically find the solution to the model (Equations ([^-([^) through simulations while vary¬ 
ing one or more parameters and keeping the others constant at the values given in Table [Tj 
We then calculate the average time spent crop-raiding once the solution has reached a steady 
state or cyclic solution. Cyclic solutions occur when baboons, ungulates and apex predators 
co-exist and, thus, to calculate the average time spent crop-raiding over the cyclic solution 
the simulation is run for 500 years and the average taken over the last 100 years. Whilst the 
long time scale is unrealistic for conservation strategies it is required to allow the transient 
dynamics to fade. Further, primates and apex predators are long-lived animals which leads to 
long period cycles so a long time scale is necessary to calculate the mean proportion of time 
spent crop-raiding. 

We use two different measures to determine crop-raiding effects, both 1 — t{B) and (1 — 
t[B))B. While only slightly different, they provide two ways to view crop-raiding dynamics - 
from the baboon perspective and from the human perspective. 1 — t{B) is the proportion of 
time a single baboon will spend crop-raiding, hence it informs us of the level of enticement for 
baboons to crop-raid. (1 —t{B))B informs us of the impact the whole population of baboons, 
crop-raiding at that level, will have on humans. It can be thought of as the number of baboons 
that choose to crop-raid 100 % of their time, while the rest of the baboons do not crop-raid at 
all. This perhaps reflects reality in that some troops of baboons do not crop-raid at all and 
other troops may solely subsist by crop-raiding. Hereafter, t{B) will be shortened to r. 


3 Results 

We will consider whether apex predator or human control is more efficient at reducing the 
proportion of time baboons spend crop-raiding, (1 — r), by focussing on four factors and their 
interactions: ( 1 ) the population of apex predators, determined by the level of human hunting 
(i7p); (2) the nutritional quality of the crops (a); (3) baboon shyness ( 7 ); and (4) human 
hunting of baboons (i7p). 


























3.1 Effect of apex predator removal on baboons and crop-raiding 
levels 

We first consider the potential for mesopredator release when apex predators are lost throngh 
overhnnting, bnt hnnting of baboons does not occnr. In Figure |^A,B,C) we show the dynamics 
of the system as hunting of apex predators is increased over time, eventually decreasing the 
apex predators to the point of extinction. This shows the striking increases in both baboon 
and ungulate populations that occur due to the loss of the apex predators. For example, the 
baboon population starts at an average of 200 baboons but increases to over 1000 as apex 
predators are driven to extinction. 

Removing apex predators, and the resulting increase in the baboon population size, has a 
knock-on effect for crop-raiding. In Figure]^ we show the proportion of time that baboons 
choose to crop-raid as the apex predator population is reduced. Due to the crowding parameter 
( 7 ), the baboons only crop-raid when they have higher population abundances, which results 
in baboons spending nearly half their time crop-raiding when apex predators are absent. 
However, the change in the effect of crop-raiding in terms of the number of baboons leaving 
the park, (1 — t)B, is even more dramatic, growing tenfold due to the loss of apex predators 
(Figure]^). When human hunting of baboons is present the graphs show a similar shape and 
pattern, but increases in population abundance of baboons and the proportion of time spent 
crop-raiding are suppressed to lower levels (see Appendix B). 

3.2 Effectiveness of human efforts to control crop-raiding 

3.2.1 Hunting of mesopredators and crop quality 

The dynamics shown in Figure assume that human hunting of baboons to discourage crop¬ 
raiding is not occurring. We now consider whether humans are able to replicate the effect of 
apex predators by hunting baboons during crop-raiding events. We examine human imposed 
baboon control both in the presence and absence of apex predators to see how much hunting 
of baboons is required to keep the level of crop-raiding low, and how it depends on apex 
predators. However, the energetic value of crops is another aspect affecting the propensity 
of the baboons to crop-raid. Higher quality crops will increase the growth rate of baboons 
more than lower quality crops, increasing the attractiveness of crop-raiding, and increasing 
the probability that baboons will raid even if they are generally shy. To understand in detail 
the amount of human hunting of baboons which is required to discourage the crop-raiding 
events, we consider how different levels of crop attractiveness affects the number of baboons 
crop-raiding. In Figure we show the difference in the number of baboons crop raiding, 
{1 — t)B, between two scenarios: apex predators present and no human hunting of baboons 
occurs versus apex predators absent and human hunting is used to control crop-raiding. We 
calculate this difference for a range of values of crop quality, a. We are able to directly see 
the change in effectiveness of the two different control strategies working separately. Positive 
values in Figure [^indicate that the number of baboons crop-raiding increases by this amount 
after apex predators are extirpated and human hunting occurs instead. For example, if we take 
the lowest black line, a = 0.5, and trace along this as Hb increases, it shows that replacing 
apex predator control with human hunting at the rate Hb will lead to increases in the number 
of baboons crop-raiding until Hb = 0 . 1 . 
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Figure 2: Gradual increase of hunting on the apex predators {Hp) to show the potential of 
mesopredator release with loss of apex predators. In A, B and C the dynamics of baboons, 
ungulates and apex predators over time, respectively. In D the proportion of time baboons 
choose to crop-raid, (1 — r), as well as the potential damage caused when they do so in E, 
measured as the number of baboons raiding, (1 — t)B, are plotted against time. Parameter 
values are as in Table[T]apart from Hp which varies from 0 to 0.3 every 250 years in increments 
of 0.06, as indicated on the top axis. The long time scale is used to show the cyclic dynamics 
at each level of hunting in more detail. 


Figure shows that virtually all values of a and Hp lead to increases in the number of 
baboons crop-raiding when subject to human rather than apex predator control, with increases 
in a more than doubling the number of baboons crop-raiding at low hunting levels. Even for 
high levels of hunting of baboons (e.g. Hp = 0.25), when crops are more nutritional {a > 1.4), 
there will always be an increase in number of baboons choosing to crop-raid 100% of the 
time when subject to human hunting rather than apex predator control. The horizontal lines 
for lower a occur because the baboons spend all their time in the park once human hunting 
is introduced (and predators are absent) as human hunting and the lack of good resources 
from crops make crop-raiding no longer worthwhile. Thus, if the crops are of poor nutritional 
value for baboons and human hunting is high, baboons will completely stop crop-raiding. In 
this case, control by apex predators leads to minimally more crop-raiding (negative values in 
Figure]^ than when apex predators are lost and human hunting is used as a control. However, 
both high levels of human hunting and poor nutritional crops (a < 1 indicates that the crops 
are less nutritional than the resources of the park) must be present for this to occur. Yet 
the baboons are predominantly still choosing to crop-raid unless hunting reaches high values; 
for example, for a = 1, Hp needs to be above approximately 0.19 to deter more crop-raiding 
under hunting than predation. This highlights the strength of apex predator over human 
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Figure 3: Varying the nutrition of the crops (a) can lead to differences in the impact of the 
crop-raiding. The increase in number of baboons crop-raiding, (1 — t)B, from when they are 
controlled purely by apex predator presence {Hp = 0, Hp = 0) to purely human hunting 
{Hp = 0.3, apex predators are extirpated) is shown against different levels of human hunting 
{Hp) for different values of a, the crop-induced increase in the growth rate, a increases from 
0.5 to 1.8 in increments of 0.1, as indicated by the arrow. All other parameters are as in Table 
The blue line indicates when there is no net increase or decrease in baboons crop-raiding. 


control in its effectiveness of reducing crop-raiding. 

3.2.2 Mesopredator shyness and crop-raiding propensity 

Shyness of the mesopredator will also affect their crop-raiding propensity. We analyze the 
potential for baboon shyness to affect crop-raiding time in Figure]^ both for apex predators 
present (Figure |^,C) and apex predators absent (Figure |^,D). We consider the impact 
of baboon shyness alongside two different management strategies, human hunting in Figure 
|^A,B) and planting crops of different nutritional value in Figure |^C,D). 

As expected, lower 7 , i.e. reduced shyness, and lower hunting levels induce more crop¬ 
raiding both for apex predators present and apex predators absent (Figure |4[\,B). When 
shyness is low crop-raiding will occur regardless of how much hunting pressure is applied. 
However, in the case when apex predators are absent (and populations are larger) baboons 
are more likely to crop-raid, including for higher shyness levels. The suppression of baboons 
by predators is amplihed by their shyness. Shy baboons will only leave the park when they 
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are forced to, due to high abundance and crowding, which corresponds to low apex predator 
numbers. The more gradual change in crop-raiding in the horizontal direction for low hfs in 
Figure compared to Figure highlights the fact that the increased density of baboons 
within the park requires 7 to be very high in order to effectively reduce crop-raiding. 

Predators No Predators 

0.5 
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Figure 4: The proportion of time baboons crop-raid (1 — r) depending on baboon shyness ( 7 ) 
about exiting the park. Baboon shyness is plotted against human hunting {Hb) in A, B and 
against the nutritional value of crops {a) in C,D. In A, C apex predators are present in the 
park {Hp = 0); in B,D apex predators are absent in the park {Hp = 0.3). In A,B, a = 1.49 
and in C, D, Hp = 0.1. All other parameters are as in Table 0 

In Figure |^C,D) we focus on the role of the nutritional value of crops (a) and mesopredator 
shyness ( 7 ), for a low level of human hunting (Hb = 0.1). In Figure]^, even with high 7 , and 
human hunting occurring, the absence of apex predators leads to baboons still choosing to 
crop-raid for a small proportion of time if the crops are benehcial enough. Again the presence 
of apex predators reduces the amount of time the baboons spend crop-raiding. However, due 
to human hunting of baboons, not all crops will coax baboons out of the park. 

It is possible to compare in Figure]^ the two strategies of human hunting and lower nutri¬ 
tional crops for different values of baboon shyness. When Hp = 0 there are higher levels of 
baboon crop-raiding than when a = 1 . 8 , indicating that keeping the nutritional value of crops 
high is of less detriment than reducing the amount of human hunting of baboons. This is true 
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for all values of baboon shyness. 


3.3 Mixed hunting strategy 

In the previous hgures we focused on the two extreme cases of no apex predator hunting 
by humans or a complete absence of apex predators (for example, from over-hunting). Now 
we wish to consider the effect of different mean abundances of predators due to changes 
in hunting impact. We examine the effects of intermediate levels of apex predator hunting 
without extirpation on crop-raiding in Figure]^ 



Figure 5: The effect of human hunting of baboons {Hb) and abundance of apex predators 
{Hp) on crop-raiding. In A the proportion of time spent crop-raiding (1 — r) is shown against 
Hb and Hp. In B the number of baboons crop-raiding, i.e. (1 — t)B is shown against Hp and 
Hp. All other parameters are as in Table 0 

As Hp increases beyond approximately 0.23 in Figure the level of hunting on apex 
predators no longer affects the proportion of time spent crop-raiding, rather it is changing 
due to Hb only. This is because apex predators die out for higher values of Hp. As Hp 
approaches 0.23, the amount of crop-raiding increases. A sharper increase occurs for the 
number of baboons crop-raiding - in Figure the color scale changes more quickly as Hp 
approaches 0.23. This is because both the time spent crop-raiding and the actual number 
of baboons increase due to loss of apex predator suppression, hence (1 — t)B increases at 
an even faster rate. This hgure highlights the potential to have control over apex predator 
abundance from hunting without reducing the benehts of apex predator suppression of the 
baboon population, since for Hp < 0.19 there are very low levels of crop-raiding occurring. 
This is true even when human hunting of baboons is at low levels. For further details showing 
the effects of this mixed hunting strategy on apex predator abundance, see Appendix B. 
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4 Discussion 


In this study we used a mathematical model to examine the trade-offs inherent in mesopredator 
release due to apex predator extirpation. We focused on an application to crop-raiding baboons 
in African parks, and the potential for crop-raiding to be mitigated by human control via 
hunting of baboons, predator population maintenance, or crop choice. 

Our results show that the presence and abundance of apex predators is one of the most 
effective strategies for reducing time spent crop-raiding by baboons as the apex predators 
suppress baboon abundance. Even when the crops are of high nutritional value, baboons 
exhibit little shyness, and no human hunting of baboons occurs, the proportion of time baboons 
crop-raid remains low as long as apex predators are present in the park to control population 
numbers. We also predict that it is possible to control both apex predator and baboon 
populations at the same time, while maintaining low levels of crop-raiding. As long as hunting 
of apex predators is kept below a certain level, apex predator numbers can be controlled 
and little hunting of baboons is necessary to keep crop-raiding at low levels. People living 
near parks may be unhappy with a high abundance of apex predators present in the park 
even if they control baboon crop-raiding, or park managers may wish to keep apex predator 
abundance under control as part of conservation strategies. Thus, a successful approach could 
be to combine the management strategies of hunting of apex predators and of baboons. 

Human hunting of baboons on its own is not able to replicate the suppression of baboons 
that occurs under predation, and so is not typically able to significantly lower the time spent 
crop-raiding in comparison. This is especially true if crops have high nutritional value and so 
induce the baboons to crop-raid. Figure [^highlights the potential for reducing the number of 
baboons crop-raiding by lowering a, the nutritional value of the crops; it is possible to halve 
the number of baboons crop-raiding when a is reduced. One potential management strategy 
could be to replace crops such as maize by less nutritious crops, such as tea. 

Higher baboon shyness leads to a reduction in the time spent crop-raiding by baboons. 
Unfortunately, baboon shyness is not an aspect of the system which the park managers or 
humans living near the park can control. However, as the shyness parameter modulates the 
amount of risk the baboons are willing to take while crop-raiding, it can also be interpreted 
as a measure of how far the baboons are willing to travel to reach the crops. If the baboons 
are very shy they would, for instance, be less willing to travel a long distance to crop-raid. 
This provides an additional type of crop control - planting attractive crops further from park 
edges. Thus, we can interpret Figure [^ with this new framework, allowing us to consider how 
three of the potential management strategies interact - less nutritious crops, planting further 
from the boundary, and human hunting. For example, if the crops are of high nutritional 
value and human hunting is occurring at medium levels, planting further away from the park 
boundary could successfully counteract the high proportion of time spent crop-raiding due 
to loss of apex predators. The requirement for both human hunting of baboons and planting 
further from the boundary in this example emphasizes the need to utilize multiple management 
strategies to control crop-raiding when apex predators are absent, which is not necessary 
when apex predators are present. At least two out of the three management strategies of less 
nutritional crops, planting further from the boundary, and human hunting must occur in order 
for management of crop-raiding to be successful if apex predators have been extirpated. 

This paper highlights the strong positive effects of top-down control by apex predators 
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on mesopredators. Mesopredator release can have wide-ranging deleterious effects. In our 
example it leads to increasing levels of crop-raiding. In other systems, trophic downgrading, 
the removal of apex predators from an ecosystem, can have far-reaching effects on species 


further down the food chain such as extinction and loss of habitat (Estes et al., 2013, 2011), 
and hence it is of great concern. These indirect effects of apex predators may not be realized 
until after the apex predators have been extirpated. This work substantiates the principle that 


all trophic levels needs to be maintained to promote healthy and balanced ecosystems (Estes 


et al. 

2011 

Ripple et al. 

2014 


Further, doing so can have knock-on benefits for humans, and 
can compliment other management strategies. The model presented here, and similar models, 
can be valuable tools in evaluating the potential impacts of suites of management strategies 
and therefore inform approaches for managing protected areas. 
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A Demography Calculations 


Using a life-table analysis, we calculated the maximum demographic potential of baboons in 
an unrestricted resource environment, using parameters from the literature. Olive baboon 
{Papio anubis) females can live 25-30 years. It is suggested that baboon females reach sexual 


maturity at 4-5 years of age and begin reproduction 1-2 years later (Altmann et al.. 1988) 


with an interbirth interval of 1-2 years (Smuts & Nicolson, 1989) although age of maturity and 


reproduction in the wild has been shown to vary with demographic and environmental factors 

They have been recorded as reproducing 


Beehner et al. 

2006 

Charpentier et al. 

2008 


as early as 5 and continuing until over 20 years of age (Smuts & Nicolson, 1989), however, 
fecundity falls off quite dramatically after this point. 


We tabulated data from Strum & Western (1982) and Altmann et al. (1977) to estimate 


fecundity and A, the rate of population growth, in a wild population to ascertain how this 
compared with unrestricted life-history calculated measures. 

We set up our life table with 30 yearly age classes (x), with a constant mortality effect of 
1/30 (1/maximum lifespan). We seeded this population with a large number (1000), so that 
rates were not influenced by small number problems. We set the age of first reproduction 
at the minimum of 5 years and termed at 25 years, to maximize reproductive lifespan. We 
then established fecundity by assuming minimum interbirth interval (1 year), with 1 successful 
offspring per event (but we chose not to include twinning), and a sex ratio skewed fully female. 
This yields a maximum A of 1.204, or a growth rate of 20.4% annually. 

There are many ways to explore the range of differences our assumptions impose on this 
estimate: for example, assuming parity sex ratio at birth, we find this drops to 14% (A = 1.14), 
because we have effectively halved female fertility. 
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Using the data from Strum & Western (1982) and Altmann et al. ( 1977[ ) we calculated 
A based on total population counts, and found a geometric mean estimate (for exponential 
growth) of 1.025, or 2.5% annual growth. In terms of female fecundity, there was a raw 
fecundity (average annual reproduction) of 0.576. 

If we combine this fecundity with our life table, assuming an equal sex ratio at birth, 
we set nix to 0.29, reducing A to 1.10, or 10% annual growth. In order to approach the 
measured A, while maintaining assumptions of constant mortality, we must impose a > 10% 
mortality rate. Clearly the life-table assumption of constant mortality and fecundity rates 
during reproductive life-span are inappropriate to real baboon demography, but this life history 
modeling exercise shows the capacity of baboons to respond to resources demographically, as 
suggested by Bercovitch & Strum (1993), in absence of any restrictions. 


B Further Analysis 

We expand on our results to show more details of how the four main parameters of interest 
affect the proportion of time that the baboons choose to crop-raid. This includes further 
information on Figure 5 in the main text and a sensitivity analysis of changes in our parameters 
on the proportion of time crop-raiding. 

B.l Effect of apex predator removal on baboons and crop-raiding 
levels 

In Figure 2 in the main text, we show how the loss of apex predators leads to a higher 
abundance of both baboons and ungulates, and this results in increases in the proportion of 
time that baboons crop-raid. Notably, there is a signihcant rise in the number of baboons 
that crop-raid, measured as (1 — t)B. In Figure we investigate the loss of apex predators 
on baboon abundance and proportion of time spent crop-raiding using the same technique as 
Figure 2 in the main text, but now with human hunting of baboons occurring. 

From Figure]^ we see that baboon abundance and proportion of time spent crop-raiding 
both increase as the apex predators are hunted to extinction. However, human hunting of 
baboons is now present which suppresses the baboon population and hence also the proportion 
of time they choose to crop-raid. Human hunting is successful at reducing the proportion of 
time crop-raiding because it is at a very high level, Hb = 0.25, which shows the potential 
for human hunting to be reasonably effective. However, even with hunting at its maximum 
level, it does not keep the proportion of time crop-raiding at the same low level as when apex 
predators are present. 

B.2 Mixed Hunting Strategy Expanded 

In Figure 5 in the main text, the proportion of time baboons crop-raid is shown against the 
combined strategies of hunting apex predators and baboons. This shows that it is possible to 
hunt apex predators up to = 0.23 and still gain the benehcial apex predator suppression 
of baboons. However, it is useful to understand what effect this hunting of apex predators has 
on apex predator abundance. Thus for each value of Hb and Hp in Figure 5 in the main text. 
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Figure 6: Gradual increase of hunting on the apex predators (Ffp) to show mesopredator 
release with loss of apex predators, alongside human hunting of baboons (Ffs). In A, B and 
C the dynamics of baboons, ungulates and apex predators over time, respectively. In D the 
proportion of time baboons choose to crop-raid, (1 —r), as well as the potential damage caused 
when they do so in E, measured as the number of baboons raiding — are plotted against 
time. Parameter values are as in Table 1 in the main text apart from Hb = 0.25 and Hp 
which varies from 0 to 0.3 every 250 years in increments of 0.06, as indicated on the top axis. 
The long time scale is used to show the cyclic dynamics at each level of hunting in more detail. 


we calculate mean apex predator abundance over the last 100 years of the 500 year simulation 
(to allow the transient dynamics to fade), as shown in Figure]^ 

From Figure it is clear that Hb does not have a signihcant effect on the number of apex 
predators as the color does not change in the vertical direction, apart from a slight curve in 
the middle region. This is because when there is high apex predator abundance baboons do 
not crop-raid. Hence changes in iJp do not impact greatly on baboon abundance as hunting 
of baboons only occurs when they crop-raid. Therefore, if baboon numbers are not affected, 
there is no feedback to changes in apex predator abundance. However, for mid-range values 
of i/p, higher hunting of baboons will allow apex predators to survive a slightly higher level 
of hunting before extirpation. Figure shows that lower levels of apex predator hunting will 
reduce the abundance of apex predators, although there will only be a signihcant decrease 
for Hp > 0.16. From Figure 5 in the main text, we know that as long as Hp < 0.23, it is 
possible to keep the proportion of time baboons crop-raid at low levels. Thus, there is a range 
of values for Hp such that apex predator abundance can be reduced whilst at the same time 
maintaining the benehts of baboon suppression. 
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Figure 7: The effect of human hunting of baboons {Hb) and apex predators {Hp) on the mean 
abundance of apex predators. All other parameters are as in Table 1 in the main text. 


B.3 Sensitivity Analysis 

In the main text, we focus on comparing one or two factors together that lead to changes in the 
proportion of time that baboons crop-raid. In Figurej^ we wish to compare the 4 main factors, 
and the two-way interactions between the factors. To do this, we followed the procedure in 


Coutts & Yokomizo (2014) by using a generalized linear model to run a global sensitivity 


analysis of the proportion of time spent crop-raiding to changes in each of the four parameters 
of interest {Hb, Hp, a and 7 ). Briefly, the method involves running 3000 simulations of the 
model while each of the four parameters can take any value within a range. We used the same 
ranges as in the main text, i.e. 0 < Hp < 0.25, 0 < Hp < 0.3, 0.5 < a < 1.8 and 0 < 7 < 1.2, 
and used Latin hypercube sampling method in order to get the most representative sets of 
parameters for the 3000 simulations. For each simulation we calculated the average proportion 
of time spent crop-raiding, 1 — r, over the last 200 out of 1000 years. After standardizing the 
parameter inputs so that they can be compared more easily, we then ran a generalized linear 
model in R using the data. As the data is in terms of proportions, we use the binomial family 
which results in the logit link function. 

From Figure we can see that a one standard deviation increase in parameters a and Hp 
leads to increases in the proportion of time spent crop-raiding, while a one standard deviation 
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Figure 8 : A global sensitivity analysis showing the effect of one standard deviation change in 
each of the 4 parameters (a, 7 , Hb and Hp) on the proportion of time spent crop-raiding. 
The black dot shows the value of the slopes and the interactions between slopes, with the grey 
bars indicating 95% confidence intervals. 


increase in Hp and 7 would lead to decreases in time spent crop-raiding. These results are in 
agreement with Figures 2-5 in the main text. If crops are more nutritional (higher a), baboons 
are encouraged out to crop-raid, as well as when their population is not being suppressed by 
apex predators (higher hunting of apex predators, Hp). Conversely, human hunting of baboons 
{Hb) and higher shyness of baboons ( 7 ) will discourage them from venturing out of the park. 
These effects are all significant; however, the only two-way interaction significantly different 
from zero is that of Hp and a. Figure [^indicates that both a and 7 are the most important to 
change in order to reduce the proportion of time spent crop-raiding. Whilst it is not possible 
to change the shyness of baboons, planting crops further from the boundary can act as a proxy 
to increase the risks a baboon would face in order to crop-raid. Thus, planting less nutritional 
crops and planting them further from the boundary could be a very successful method to 
reduce crop-raiding. However, the benefits of this would have to be outweighed against the 
costs of different crop strategies and loss of usable farmland. The hunting of apex predators, 
Hp is not shown to have as influential an impact as a in Figure This is likely caused by the 
fact that Hp has a threshold effect. Once Hp is large enough the apex predators go extinct, 
thus further increases in Hp will not affect the dynamics. Therefore comparison between 
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the sensitivity of the model to changes in Hp and the other parameters should be inferred 
with caution. From the main text, we know that the presence of apex predators significantly 
reduces the proportion of time spent crop-raiding. 
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